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Volume element structure and roton-maxon-phonon excitations in superfluid helium 

beyond the Gross-Pitaevskii approximation 
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We propose a theory which deals with the structure and interactions of volume elements in 
liquid helium II. The approach consists of two nested models linked via parametric space. The 
short-wavelength part describes the interior structure of the fluid element using a non-perturbative 
approach based on the logarithmic wave equation; it suggests the Gaussian-like behaviour of the 
element's interior density and interparticle interaction potential. The long- wavelength part is the 
quantum many-body theory of such elements which deals with their dynamics and interactions. 
Our approach leads to a unified description of the phonon, maxon and roton excitations, and has 
noteworthy agreement with experiment: with one essential parameter to fit we reproduce at high 
accuracy not only the roton minimum but also the neighboring local maximum as well as the sound 
velocity and structure factor. 
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INTRODUCTION 



The microscopical structure of quantum liquids is dif- 
ferent from that of classical ones. The latter are discrete 
fluids consisting of localized atoms or molecules which 
thermal de Broglie length is smaller than the average 
atomic or molecular separation. For instance, the ther- 
mal de Broglie wavelength of a water molecule at room 
temperature is less than the radius of the hydrogen atom. 
Therefore, classical liquids can be treated as continuous 
media only on the length scales which are larger than 
atomic or molecular separation. In quantum liquids it is 
other way around: atoms become delocalized in space be- 
cause their de Broglie lengths are always larger than the 
inter-atomic distance, and overlapping with each other. 
As a result, instead of discrete atoms, quantum liquids 
must be described in terms of the fluid volume elements 
and elementary excitations such that non-locality and 
continuity are preserved down to very short scales 0- 
3]. Historically, it was Landau [|| who developed the 
theory of superfluidity based on the mechanics of contin- 
uous media which qualitatively agreed with experiment 

However, some questions remain open in this regard - 
in particular, those about the interplay between original 
degrees of freedom ( 4 He atoms) and emergent collective 
ones - such as the fluid volume elements (also called the 
fluid particles in the Lagrangian description Q). Indeed, 
the fluid-dynamical description of liquid helium presumes 
that usage of the volume element's notion must be phys- 
ically justified before, since it is the volume element, not 
the atom itself, which is supposed to be the most ele- 
mentary object of the fluid approach. Thus, the under- 
standing and describing of this notion in the quantum 
fluid theory is more important than in the classical one 
because for classical fluids the continuity is just the long- 
wavelength approximation whereas for quantum ones it 
is an empirically established fact. Therefore, searches for 



the most suitable collective degrees of freedom in super- 
fluids continue, and the issue of how the formation and 
stability of the volume element of liquid helium as a con- 
tinuous medium can be explained from first principles 
remains a subject of active studies nowadays [3. [To|. 

The Bose-Einstein condensation (BEC) is another as- 
pect of the superfluidity phenomenon which must be 
taken into account [Til ]. Being predicted almost a cen- 
tury ago, it was used by London and Tisza |"T5 to explain 
the superfluidity of liquid helium discovered by Kapitsa 
and Allen [l3|]. The existence of BEC in the superfluid 
phase has been confirmed by inelastic neutron scattering 
|l4l [lH. thus, the influence of BEC upon the properties 
of liquid helium is under current study as well [ljH18[ . 

Here we describe the structure of the volume element 
of superfluid helium as well as the phcnomcnological con- 
sequences which follow. Among the derived effects is the 
famous Landau roton spectrum for which full theoretical 
explanation is still pending despite the numerous efforts 
made towards its understanding [1MG1I3. As a mat- 
ter of fact, even usage of the term "roton" itself is often a 
purely historical one since the explanation of the Landau 
spectrum in helium does not necessarily have to be based 
on rotating degrees of freedom of helium atoms 0, . 
Indeed, here we show that it is possible to formulate a 
theory without their explicit involvement - yet obtain a 
remarkable agreement with experiment. 

Throughout the paper we neglect temperature effects 
which is a standard assumption for fundamental mod- 
els since below the lambda point the thermal effects can 
be considered as corrections (the superfluid helium has 
almost perfect heat conduction). 

In section [2] we formulate a general theory of helium 
II, in sections [3] and [4] we derive the observable quanti- 
ties and compare them with available experimental data. 
Conclusions are made in Sec. [SJ 
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2. COLLECTIVE VARIABLE THEORY 

We divide a theory into two nested parts. The shorter- 
wavelength part describes the interior structure of the 
volume element of helium superfluid using the non- 
linear logarithmic quantum wave equation. The longer- 
wavelength part is the quantum many-body theory of 
these elements, treated as new collective degrees of free- 
dom. While one assumes them to be effectively point- 
like objects in the long-wavelength approximation, their 
spatial extent and internal structure are taken into ac- 
count by virtue of the nonlocal interaction term. At that, 
these two models are not independent: in the shorter- 
wavelength part we derive behavior and specific values 
of the parameters for the longer-scale one, in accordance 
with the "Russian-doll" picture of nesting of scales. 

As a starting point, we assume that strongly-correlated 
helium atoms can form a bound state characterized by 
a single macroscopical wavefunction. The wave equation 
describing such object cannot be of the Gross-Pitaevskii 
(cubic Schrodinger) type for at least two reasons. 

The first one is that the GP approach 21] is a per- 
turbative one which takes into account only two-body 
interactions and neglects anomalous contributions to self- 
energy which is a good approximation for dilute system 
like cold gases [22], but unlikely to suffice for liquids. In- 
deed, according to aforesaid the atoms in quantum liquid 
are delocalized and thus nothing prevents them from get- 
ting involved into multiple-body interactions where the 
multiplicity can vary from two to the total number of 
atoms in the system. An example of multi-body (three 
and more) interactions being very important for forming 
bound states of bosons at low temperatures is the Efimov 
state I2 3| . | 2 4ll which has been experimentally observed in 
helium|25l|7a recent review can be found in [26j . 

The second obstacle is that the ground-state wavefunc- 
tion of the free-space GP BEC model does not describe a 
localized object. Instead, the free GP condensate tends 
to occupy all available volume - as such one needs to ap- 
ply an external potential trap to confine the condensate 
and stabilize the system. It is difficult to imagine that 
upon the transition into a superfluid phase the free he- 
lium suddenly becomes surrounded, both in parts and 
as a whole, by a hypothetical external potential, not to 
mention that the coherent appearance, stability and syn- 
chronization of such domains across the bulk would be 
impeded by significant volatility of the liquid. 

There exists, however, another candidate where these 
problems simply do not occur in first place. This is non- 
linear Bose liquid defined by virtue of the logarithmic 
Schrodinger equation: 



h 2 -. 
2m 



^ 1 ln(5 3 |*| 



* = 0, (1) 



where — ^(x, t) is the wavefunction of condensate nor- 
malized to the number of particles N - such that parti- 
cle density is determined as n = \^\ 2 , m is the mass of 
the constituent particle (helium-4 atom in our case, i.e., 



m = itihc ~ 6.64 x 10~ 24 g), and f3 and a are constant pa- 
rameters of interaction. The equation alone received at- 
tention a while ago as the simplest U(l)-symmetric wave 
equation (apart from the conventional Schrodinger one) 
which satisfies the dilatation covariance and separability 
properties [13, HH but for the purposes of theory of quan- 
tum liquids it has been applied only recently [29| . While 
the second-quantized Hamiltonian which might lead to 
such equation is not known yet, the equation itself has 
some interesting properties and thus can be used as a 
starting point. The equation can be also derived from 
the field-theoretical action with the Lagrangian 

C = y (<W* - V*d t V) + |^|W| 2 + XM|*| 2 ), (2) 



where the potential is defined as 

Vf}{n) = -p- lr n [In (n« 3 ) - l] 



(3) 



The latter opens down and has local non-zero maxima at 
"oxt = 5~ 3 for positive j3, see Fig. [T] Despite the poten- 
tial not being bounded from below as a function of \^/\, 
no density divergences arise since the wavefunction can- 
not take arbitrarily large values, due to the normalization 
constraint. 
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FIG. 1: Field potential for the logarithmic BEC (in units 
/3d 3 ) versus square root of particle density (in units a~ 3 ^ 2 ). 

It turns out that the logarithmic Bose liquid has a 
number of features suitable for our objectives: it im- 
plicates not only binary but also multiple-body interac- 
tions (when more than two bodies can scatter simultane- 
ously), and its ground state is the so-called gausson [28j 
- a spherically-symmetric object which is localized and 
stable even in absence of a trapping potential, with the 
interior density obeying the Gaussian law 



n{x) = n(0)e" (r/a)2 , 



(4) 



where n(0) = A/(?r 3 / 2 a 3 ) is the central particle density, 
r = |aJ|, and a — h\J j3 /2m is the characteristic radius. 
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In principle, for practical purposes the relation a ~ a 
can be assumed throughout this paper, which essentially 
means that the parameters of the short-scale model are 
not entirely independent but must be bound by (at least) 
one physical constraint, rs p 2 /2m, where p — h/a is 
the de Broglie momentum corresponding to the length 
scale a. 

It should be emphasized that the object described by 
((4]) is different from the classical droplet since it does 
not have border in a classical sense, therefore, its stabil- 
ity is supported not by surface tension but by nonlinear 
quantum effects in the bulk (29|. 

The long- wavelength part can be formulated as follows. 
As long as we have established that the nucleation and 
stability of volume elements of the Gaussian type is justi- 
fied on quantum-mechanical grounds, it is natural to as- 
sume that below the critical temperature the atoms tend 
to form the fluid elements of this kind. Therefore, new 
collective degrees of freedom arise, so a theory must be 
formulated in terms of the volume elements rather than 
of the helium atoms themselves [IcJ. Besides, as long 
as the bosonic liquid can no longer be assumed homo- 
geneous one cannot use all the results of section 4 from 
29]. Thus, we use the following Hamiltonian 



ip{x) 



H-, 



int. . 



(5) 



where M is the mass of the volume element and ip is the 
corresponding field operator. The interaction is defined 
via the nonlocal term 



Hint — 5 



d 3 xd z x'^{x)^{x')U{\x-x'\)^{x')^{x), 



where U(r) is the energy of interaction between volume 
elements. This energy can be estimated in the following 
way. Using (jj) one can derive that the Gaussian volume 
element of size r ~ a stores an amount of internal bulk 
mass-energy 



f 



e(r) oc / n(r')r' 2 dr' oc -(r-r ) e^ r/a) [1 + O (r - a)} , 
Jo a 

where the value r = a [1/2 + l/(ey/ired{l))] « 0.75a 
refers to the point where the dominant term of e(r) 
changes sign. Since each element has been shown to 
be stable with respect to small perturbations, it tries to 
maintain its size and mass when interacting with immedi- 
ate environment, therefore, to alter these values one has 
to supply the amount of energy which is proportional to 
e(r). In absence of strong external fields this energy can 
come only via interaction with other volume elements 
hence we can conclude that U(\x — x'\) oc e(\x — x'\). 
Thus wc can introduce the proportionality factor £/o(i?), 
assume it to be a constant in a leading-order approxima- 
tion, and estimate the interparticle potential as: 



U(r) = ^( r -r )e- (r/a)2 , 



(6) 



up to the terms of order O (r — a) which are assumed 
to be small. These terms can be safely omitted unless 
interaction deforms the elements so strongly that their 
interior structure cannot be neglected anymore; but then 
the element's spherical symmetry becomes deformed as 
well. The quantity U = -aU(0)/r -1.34 U(0) be- 
comes the free parameter of the long-wavelength part of 
the theory. If Uq is positive then the critical radius ro de- 
termines the inter-element separation below which a pair 
of neighboring volume elements becomes unstable against 
coalescence. In principle, any possible effects of the de- 
formations of elements can be accounted for by further 
upgrading the constant Uq to a correspondingly derived 
function Uo(x) but for the purposes of our current study 
the approximation Uq(x) ~ Uq will do the job, as will be 
shown below. 



3. OBSERVABLES 

In this section we derive the analytical expressions for 
the main observables of the theory: energy of excitations, 
structure factor and speed of sound. 



3.1. Energy of excitations 

To derive the energy spectrum of excitations one can 
use three alternative approaches. The first one is based 
on the perturbation theory [la] . second is the standard 
Bogoliubov transform [30l ]. and third is about analyzing 
small perturbations of the equation of motion for the field 
operator - as shown in [l8| by the example of the semi- 
transparent sphere model. In our case one can check 
that all three approaches yield the same spectrum in the 
leading approximation. The Bogoliubov's method is the 
most straightforward and simple one, hence it will be 
used from now on. After some algebra we arrive at the 
following ./V-body Hamiltonian operator 



H « ^nlU V + E n ala r , 



(7) 



where ap and a p are the creation and annihilation opera- 
tors of the quasi-particle with momentum p, hq — N/V — 
pHe/M is the background particle density of the liquid of 
N elements occupying the volume V, pHo ~ 0.145 g/cm 3 
is the liquid helium-4 density, U p = J U (x) e l ^' x / h d 3 x is 
the Fourier transform of the interaction potential ©, and 
E p is the quasi-particle's energy counted from a ground 
state. The energy obeys the following dispersion relation: 



E„ = -i— Jl 
p 2M V 



4n MU p 



with the Fourier transform being calculated as 

,3; 



U p = ira^Uo (l - V^fi 



ke -(ak/2Y 



(8) 



(9) 



where fk = ro/a+ \_\ak — (ak) erfi(afc/2) and k = p/h 
is the wave vector. 
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FIG. 2: Energy of quasi-particles E p versus momentum (in 
units of E a and p a , respectively) for different values of ito- 
Below the value 11.43 the roton minimum disappears, and 
above 25.98 energy becomes complex-valued. 

Introducing the volume element's de Broglie momen- 
tum scale p a = 2%/a and corresponding energy scale 
E a = p a /2M, one can re- write the dispersion relation 
(|SJ) in the dimensionless form which shows that the dy- 
namics of the quasi-particle depends on just one essential 
(non-scale) parameter 

w = ZimQcfMUa/h 2 . (10) 

Numerical analysis shows that admissible values of this 
parameter lie between 11.43 and 25.98 where the local 
(so-called roton) minimum appears, see Fig. [2l In its 
vicinity the dispersion relation can be written in the Lan- 
dau form, E p w A+ ^~ Po ^ , where A is the roton energy 
gap, po is the minimum value of momentum, fj, is the 
effective roton mass. The corresponding profiles for dif- 
ferent values of uq are given in Fig. [3] from where one 
can see that the roton mass diverges at the lower bound 
(the local minimum becomes a saddle point) whereas the 
energy becomes complex-valued above the upper bound. 

We can also plot the quasi-particle's velocity v = vl 
where v — dE p /dp and I — p/p. From Fig. |4]one can 
see that within the above-mentioned range of the param- 
eter uq the liquid indeed exhibits the superfluidity fea- 
ture - an interval of momentum and energy over which 
the effective mass p/v turns negative always exists, thus 
indicating that the production of dissipative excitations 
is suppressed. As a matter of fact, this is just a restate- 
ment of the original Landau idea of introducing the roton 
minimum into the excitation spectrum to explain the su- 
perfluidity - since a negative p/v range always appears if 
a positive-definite function has the local minimum pre- 
ceded by a local maximum. One can also notice that at 
small momenta (phonon regime) the velocity varies very 
slowly as a function of frequency uj = E p /H, as if there is 
no dispersion at all. 
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FIG. 3: Local minimum's values of momentum po/Pa (dashed 
curve), roton energy gap A/E a (solid) and mass n/M (dotted) 
versus the interaction coupling constant taken between mo ~ 
11.43 and 25.98. 

3.2. Speed of sound 

By sound here we understand the conventional acous- 
tic oscillations, the heat transfer and second sound are 
not considered since the temperature effects are neglected 
throughout the paper, as mentioned in the introduction. 
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FIG. 4: Quasi-particle's velocity v (in units of w| p =o) versus 
momentum (in units of po, dashed curve) and energy (in units 
of A, solid curve). It is plotted at uo ~ 20.6 (then «| p= o ~ 273 
m/s) but such behavior holds qualitatively for uq between 
11.43 and 25.98. 

Thus, our observable of interest is a mean value of the 
sound velocity (v s ) which can be estimated as follows. 
If size and interior structure of volume elements are ne- 
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glected then the speed of sound would be given just by 
the standard formula 12 211 



,(°) = 



uqUq/M = lim v 



(11) 



However, in general one should take into account that 
sound propagation is affected by the volume element's 
"interior" which is described by the logarithmic Bose liq- 
uid in a ground state, and speed of sound inside the latter 



is known to be 



1 



(12) 



in the leading approximation [29j. The mean speed of 
sound thus can be estimated as an averaged sum of the 
background and excitation contributions taken with the 
weights controlled by the interior density of a volume 
element ((3]), 



\c s + (l-x) v 



(0) 



dr 



,(°) 



(13) 



where x — (1/2) n(x)/n(0) = ^e ( r / a )" is the weight 
factor, |7ra 3 = 1/no = M/pB. e is the effective packing 
volume of an element. The numeric coefficient inside the 
weight factor has been chosen in such a way as to satisfy 
the following limit case properties: if a — > then the vi°^ 
term dominates (zero-size limit), and if a — > oo then both 
terms contribute equally. 



3.3. Structure factor 

According to (|4|, our volume elements are essentially 
quantum states of delocalized atoms which do not have 
a border or interface in a classical sense. Therefore, their 
presence in the liquid helium can be judged only indi- 
rectly - via the structure factor Sk which determines the 
scattering of neutrons (or x-rays, after multiplication by 
the atomic structure factor) and can be measured exper- 
imentally. This situation, however, is not something ex- 
traordinary: the microscopical objects employed in pre- 
vious models of liquid helium, such as vortex rings [19j , 
hard spheres [13, EH or stochastic clusters [|[ , were never 
claimed to be "directly seen" (whatever it could mean in 
our essentially quantum context) either. Instead, if a 
selected model's predictions, e.g. for the excitation en- 
ergy and structure factor, coincide with experimentally 
observed values then the corresponding microscopical in- 
terpretation can be adopted - at least, until the better 
one arrives. 

The computation of the structure factor is slightly 
more difficult than that of the energy of excitations: as 
shown in the previous section, the computation of E p re- 
quires, basically, the knowledge of the inter-element inter- 
action potential, such as ([S]), only. By contrast, to derive 
the structure factor one should derive the excited-state 
wavefunction of the liquid helium which can be done only 
approximately under certain simplifying assumptions. If 
one considers the "bare" helium atoms to be the fun- 



damental degrees of freedom then the structure factor 
would be given by the formula derived in (l9| . In our 
case it would yield 

S ™= P 72M Ep =(l + ^)~ ,/2 , (14) 

where all quantities have been defined in the previous sec- 
tion. However, this formula can not be blindly applied 
to our case since we use the collective degrees of free- 
dom which are not "bare" helium atoms but their bound 
states. Besides, the logarithmic non-linearity affects not 
only the energy of excitations but also the correlation 
functions themselves. Therefore, one needs to slightly re- 
vise the Feynman's derivation by taking these effects into 
account - especially if one considers that Feynman and 
Cohen's theory gives only a qualitative agreement with 
experiment. Fortunately, the computation is straight- 
forward and necessary corrections are minimal (yet im- 
portant). We begin by assuming that the total energy 
of the system of N fluid elements comes from minimiz- 
ing the integral £ = J i^"Ki\) d 3N x, subject to the fixed 
normalization integral J — J ip*ip d 3N x, where ip is the 
excited-state wavefunction. Throughout this section we 
use x N to denote the set of coordinates Xi of all the el- 
ements, and J d 3N x to represent the integral over all of 
them. The Hamiltonian of the system is given by a stan- 
dard expression, 

i 

where V is the potential energy of the system and Latin 
indices run from 1 to N. Following Feynman, for the 
excited-state wavefunction we assume the Bijl ansatz 

& 

^ = i^o = £/(£iM>, m=e its , (16) 
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with the only difference that ipo = ipo(x N ) is now the 
ground-state wavefunction of one-particle states of loga- 
rithmic quantum liquid (not plain atoms as was in the 
original Feynman derivation). Due to the separability 
property of the logarithmic Schrodinger equation, they 
satisfy its stationary many-body analogue: 

Ho^o = [NE + /T 1 In (a 3A Vo| 2 )] </>o, (17) 

where the logarithmic term has replaced the potential 
part, according to the two-scale structure of our the- 



ory (we remind also that throughout the paper the ex- 
ternal potential is neglected). With the use of the 
separability property the TV-particle equation (fT7|) can 
be analytically solved for the case of a ground state, 
cf. section [2j One thus obtains the eigenvalue Eq = 
3/3- 1 [1 + ln(vWa)] ~ 3/3" 1 (l + ± In tt) and oc 

n*(2i) ocn«p(-i^i 2 /2o 2 ), cf. nf. 

i i 

Further, with the ansatz (|16p in hand the energy per 
particle of the system can be approximately written as 



E p + E = 



h z l 

1M N 



E V^* • ViF + § \F\ 2 (e + 1 In (~a 3 P T) 



p N d 3N x 



where 



Pn 



1 



/ \F\ 2 PN <P N x 



N v— v ■ — * ■ 2 / 2 
E Fi| /a 



7r 3 / 2 a 3 



(18) 



(19) 



is the ground-state density function. By construction, excitation energy E p takes into account also the inter-element 
interaction, see Sec. 13.11 Further, using the fact that logarithm varies significantly slower than the Gaussian, we 
obtain 



E p + Eq 



M 



Eo+p-'hiia^Tm 



-2 ^E/ViF* -V t F PN d 3N x 



2M J \F\ 2 PN d? 



N,, 



which after integration reduces to 



h 2 j'Wf*(x) -Wf(x)d 3 x 

P+ ~ 2M//*(fi)/(f 2 )P(fi ~x 2 )d 3 xid*x 2 ' 



(20) 



(21) 



where P(x\ — x 2 ) is the probability of finding an element at x 2 per unit volume if one is known to be at X\, and we 
denoted 



Sw£; o-^:(3 + ln^« 3/T 1 



m If 2m, 
1 h - 1 | In 7r 

M 2\ 3M 



r 



The variation with respect to /* yields the equation 

(E p + E) J f(x 2 )P(xi ~ x 2 )d 3 x 2 + ^-V 2 /(£i) = 0, 

(22) 

and upon remembering (|16l) this gives us the final for- 
mula for the structure factor Sk = J P(x)e ikve d 3 x in a 
logarithmic theory: 

-l 



S k = 



P 



2M {E p + E) 



'1 



AnaMUr, 2ME 



lQMU p 

V 



(23) 

Thus, the only difference from the "non-logarithmic" for- 
mula (|14p is the appearance of an additive constant term 



in the denominator which is induced by the logarithmic 
nonlinearity. 



4. THEORY VERSUS EXPERIMENT 

In this section we compare the derived observables of 
our model with available experimental data. If we choose 
the value uq rs 20.6 then the established experimental 
data for energy of excitations (see, for instance ||) can 
be successfully fit. When fitting, we do not use the least- 
squares techniques but rather aim at the precise position 
of the roton minimum - so as to estimate the magnitude 
of deviations from experiment at the neighboring local 
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maximum £" max at p = p ma . x also known as the maxon 
peak. It turns out that our theory can reproduce the 
maximum with accuracy which is better than three per- 
cent, see the cumulative comparison between theory and 
experiment given in Fig[5] and Table HI The other val- 
ues of primary and secondary parameters of the theory 
computed at uq ~ 20.6 turn out to be the following: 

a ps 2.37 A, M/m ps 1.228, U /A ps 69, 

a ps 1.25 A, p a /p Q ps 0.84, E a /A ps 1.46, (24) 

p- x /k B ps 3.84 K, /T 1 /^ « 0.31, 

hence this set will be used for further comparisons with 
experiments in this paper. For instance, using (|24|) the 
velocity- related quantities from section ^. 2l are evaluated 
as c s ps 89 m/s and uj, ps 273 m/s, therefore, the mean 
value of the speed of sound (flU)) can be estimated as 

(v s ) ps 231 m/s, (25) 
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FIG. 5: Energy of quasi-particles iJj, versus wave vector. Solid 
curve is a theoretical prediction, dots denote experimental 
data [eJ. 

Finally, the structure factor (|23|) is also evaluated us- 
ing the values (IMl) . We compare the theoretical curve 
with the neutron scattering data in Fig. [5] (the neutron 
scattering is the most relevant method here since it is the 
one which is being used for measuring E p as well). One 
can see that the asymptotic behaviour of the structure 
factor, both at small and large momenta, is compatible 
with known experimental data, see also the discussion in 
the concluding section. For instance, from the scattering 
data one can not surely conclude that at small momenta 
Sk must tend to zero as a linear function of p. Indeed, 
as one can see from Q (see also the first figure from 
the Feynman- Cohen paper in [ll|), Sk is not a straight 
line near origin but rather a curve which resembles more 
of a parabola and thus agrees with the asymptotic be- 
haviour of (|23l) at p — > 0. Though, one can not a priori 



exclude that in real scattering experiments some side ef- 
fects might appear and introduce the asymptotic correc- 
tions proportional to p and thus override the non-linear 
asymptotics; however, these additional effects are not a 
subject of our study at this stage. 
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FIG. 6: Structure factor Sk versus wave vector. Solid curve is 
a theoretical prediction given by (|23[) evaluated at the values 
(|24p , triangles denote the experimental data [?J ■ For the sake 
of comparison, the result of the formula (|14|l is being plotted 
as well (dashed curve). 



5. CONCLUSION 

The superfluidity of liquid helium II is the complex 
phenomenon exhibiting a unique feature of the Landau 
spectrum of excitations which successfully explains the 
behavior of specific heat, viscosity and sound velocity. 
The microscopical theory of this phenomenon must ad- 
dress many different aspects such as the Bose-Einstein 
condensation, appearance of collective degrees of free- 
dom, non-locality and continuity, formation of the vol- 
ume element of the liquid, its structure and stability. 

Here we proposed a concise analytical theory of struc- 
ture and excitations in superfluid helium which elabo- 
rates on these aspects. It consists of two models which act 
on different length scales, but are connected via the para- 
metric space: the behaviour of quantities and the values 



TABLE I: Theory versus experiment 



Quantity 


Theory 


Experiment 


Difference, % 


A/fc s , K 


8.65 (fit) 


8.65L0.04 


0±0.4 


po/h, A" 1 


1.91 (fit) 


1.91L0.01 


0±0.5 




0.16 


0.16L0.01 


0±6 


-Emax/fcs, K 


14.3 


13.92L0.1 


2.7=f0.7 




1.08 


1.11L0.04 


2.6±3.6 


(v 3 ), m/s 


231 


237±2 


2.6±0.08 
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of parameters in the long-wavelength model are derived 
from the short- wavelength part. The latter advocates 
the appearance of the collective degrees of freedom which 
justify the possibility of the fluid-dynamical description 
and can be used for characterize the volume elements of 
the quantum liquid. The "interior" structure of a fluid 
volume element is described using the non-perturbative 
approach based on the logarithmic wave equation. It 
turns out that the interior density of the element obeys 
the Gaussian distribution. 

The long-wavelength part is the quantum many-body 
theory of the volume elements as effectively point-like 
objects - yet their spatial extent and internal structure 
are taken into account by virtue of the nonlocal inter- 
action term. The corresponding interparticle interac- 
tion potential © was not postulated but derived from 
the short-wavelength part, and it appeared to be very 
simple as compared to the empirical and semi-empirical 
inter-molecular potentials widely used in perturbative 
approaches (20|. On the other hand, more fundamen- 
tal approaches, such as those based on vortex rings fl9j . 
hard spheres [T3, [H| or stochastic clusters either have 
shown so far only the qualitative agreement with exper- 
iment, to our best knowledge, or some of observables 
in those models were not computed. In our theory the 
quantitative agreement with various experimental data 
has been achieved: with only one essential parameter in 
hand, uq, we reproduced at high accuracy (better than 
three per cent) not only the roton minimum but also the 
neighboring local (maxon) maximum. The velocity of 
sound and structure factor are also computed and found 
to be in a very good agreement with experiment. Slight 
deviations of the energy spectrum and structure factor 



curves from experimentally observed ones at large values 
of momenta p > 2h A -1 come due to the scattering parti- 
cles begin to probe the interior of fluid elements and thus 
affect the corresponding wavcfunctions. At even higher 
momenta, p > Ah A" 1 , the agreement will eventually re- 
cover since the contributions from kinetic energy come to 
predominate those from interaction. 

To conclude, the collective variables we have been 
using ultimately provide a unified description of the 
phonon, maxon and roton excitations. In future it would 
be interesting to study the applications of the mod- 
els with the non-polynomial ground-state wave equa- 
tions and Gaussian-like interparticle potentials, like ([T]) 
and ([6]), in the physical situations when superfluid be- 
comes effectively low-dimensional (cigar- or disk-shaped) 
and/or subjected to external fields or admixtures. A par- 
allel direction of research would be to take into account 
temperature effects (see the introductory section), such 
as heat transfer, second sound, etc., which is necessary 
for a number of practical applications. 
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